Galactic Phase Transition at Ec = 0.11 eF from Rotation Curves of Cored LSB 
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We analyze the a set of seventeen rotation curves of Low Surface Brightness (LSB) galaxies from 
the The HI Nearby Galaxy Survey (THINGS) with different mass models to study the core structure 
and to determine a phase transition energy scale {Ec) between hot and cold dark matter, due to 
nonperturbative effects in the Bound Dark Matter (BDM) model. Our results agree with previous 
ones implying the cored profiles are preferred over the N-body motivated cuspy Navarro-Frenk- 
White (NFW) profile. We find an average galactic core radius of Vc = 260 x 10 pc and a phase 
transition energy _Ec = 0.11 x IQ*'''''' eV, that is of the same order of magnitude as the sum of the 
neutrino masses. 



Introduction. Understanding the distribution of 
dark matter (DM) in galaxies has been a major work 
in recent times [l|, [^ . Late- type Low Surface Galaxies 
(LSB) galaxies are of special interest since it is believed 
that they arc dominated by DM. and high resolution HI , 
Ha and optical data can help to distinguish among the 
different DM profiles proposed in the literature. There 
are essentially two types of halo profiles, the ones stem- 
ming from cosmological A^-body simulations that have 
a cusp in its inner region, e.g. Navarro-Frenk- White 
(NFW) profile 0]. On the other hand, the phcnomcno- 
logical motivated cored profiles, such as the Burkcrt or 
Pseudo-Isothermal (ISO) profiles Q. Cuspy and cored 
profiles can both fit most LSB rotation curves, such as 
[a, Q , where LSB rotation curves has been analyzed with 
the most relevant profiles, showing a preference for a 
cored inner region with constant density and the stel- 
lar mass models is relevant in determining the size of the 
core radius Tc- Different systematics may play an impor- 
tant role in the observations such as noncircular motions, 
resolution of data and other issues [a, Q : as well as the 
type of galaxies involved Q . There are attempts to recon- 
cile both approaches through evolution of DM halo pro- 
files including baryonic processes [lO[ to transform cuspy 
to shallower profiles that follow the solid-body velocity 
curve {v ~ r) observed in late-type LSB. However, this 
issue is a matter of recent debate [11|. Recently one of 
us proposed a new t ype of Dark Matter, called bound 
dark matter (BDM) [l2|, motivated by particle physics 
in which DM particles are relativistic at high energy den- 
sities, pbdm > Pc above the phase transition scale energy 
scale Pc ~ E'^, i.e. we have Hot DM (HDM) with pecu- 
liar velocity v ~ c. But, for lower energy densities than 
Pc the BDM particles acquire a large mass due to non- 
perturbative physics and they behave as standard Cold 
DM (CDM) particles with v <^ c. The phase transition 
between CDM and HDM given by the scale pc = Ef 
can be determined theoretically or phenomenologically 
by consistency with cosmological or galactic data. In the 
present work we will estimate its value through the study 



of the rotation curves of LSB galaxies. 
The work here is twofold. On the one hand, we use galac- 
tic rotation curves to extract information on the possi- 
ble core nature of DM dominated galaxies, and on the 
other, we determine the fundamental parameters behind 
our DM proposal. To perform these tasks we use four 
different halo mass profiles (NFW, BDM, Burkert, ISO) 
and five stellar mass models (Min.Disk, Min.Disk+gas, 
Kroupa, diet-Salpeter, Free T^,). We present here the 
main results for the BDM and NFW profiles, and we refer 
to the most extensive work in [l3|. The standard devia- 
tion of the parameters, e.g. the core radius Tc, depends 
on the stellar mass model used but the results are con- 
sistent with each other. In some galaxies the core radius 
may be different from zero at %95 confidence level and 
consistent with Tc = for different stellar mass model. 
This clearly shows the importance of having better inner 
data and a better understranding on the mass models. 
For the analysis we use The HI Nearby Galaxy Sur- 
vey (THINGS) [IJI, which collects high resolution and 
excellent sensitivity of the velocity fields revealing ex- 
tended measurements of gas rotation velocities and cir- 
cular baryonic matter trajectories [Tj]. We will show that 
the LSB rotation curves yield a phase transition energy 
scale Ec, between HDM and CDM for our BDM profile, 
cV. This Ec is a new fundamental 
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scale for DM and can also be theoretically determined 
using gauge group dynamics. However, even though we 
propose Ec as a new fundamental constant for DM it is 
important to notice that its value does depend on the 
choice of BDM profile used and on the quality of the ob- 
servational data. The coincidence in the size of the sum 
of neutrino masses with the magnitude as Ec could open 
an interesting connection between the generation of DM 
and neutrinos masses. 

BDM ModeL Cosmological evolution of gauge groups, 
similar to QCD, have been studied to understand the 
nature of dark energy [T^ and also DM [l^. For an 
asymptotically free gauge group the strength of the fun- 
damental interaction increases with decreasing energy 



and the non-pcrturbative mechanisni generates the mass 
of bound states, as baryons in QCD. In this case, the 
mass of the bound states particles is not the sum of its 
component particles but it is due to the binding energy 
and is parameterized by Ac. The condensation or phase 
transition scale is defined as the energy where the gauge 
coupling constant g becomes strong, i.e. 5(A) 3> 1, giv- 
ing a condensation scale Ac = A; e^'^'^ /''s- , where b is 
the one-loop beta function which depends only on the 
number of fields in the gauge group (for example for 
a SUSY gauge group SU{Nc),Nf, where Nc{Nf) is the 
number of colors (flavors), we have b — 3Nc — Nf) and 
gi is the value of the coupling constant at an initial scale 
Ai. Clearly, Ac is exponentially suppressed compared 
to Ai and we can understand why Ac is much smaller 
then the initial A^, which may be identified with the 
Planck, Inflation, or Unification scale. The order of mag- 
nitude of the mass of these particles is msDM = dEc 
with d = 0(1) a proportionality constant. In QCD 
one has Aqcd — 200 MeV with the pion (proton) mass 
ttItt — 140 MeV {irib — 939 MeV), giving a propor- 
tionality constant is in the range 0.7 < d < 5, with 
bound mass much larger than the mass of the quarks 
(m„ ~ (1 - 3) MeV,md ~ (3.5 - 6) MeV). In our case 
the gauge group and elementary fields are not part of 
the standard model (SM). Our dark gauge group is as- 
sumed to interact with the SM only through gravity and 
is widely predicted by extensions of the SM, such as brane 
or string theories. Wc can relate Ec to Ac since the en- 
ergy density depends on the average energy per particle 
and the particle number density n, i.e. pc = E^ = Ac n. 
There are two natural places where one may encounter 
high energy densities for dark matter. One is at early cos- 
mological times and the second place is in galactic inner 
regions. Here we are interested in the latter case. Away 
from the center of the galaxies the energy density de- 
creases and one has pbdm < Pc- In this region the BDM 
particles are CDM. Since CDM is well parameterized by 
a cuspy NFW profile pnfw = Po/[^/?'s(l + f/rs)'^)] we 
expect BDM to have this limit away from the galactic 
center and as long as pbdm ^ Pc- However, in the inner 
region the DM energy density increases and once pbdm 
reaches pc we have a phase transition and the BDM par- 
ticles become relativistic with a dispersion velocity u ~ c 
and therefore forming a core inner region. At this stage 
the NFW profile would no longer describe the behavior 
of our BDM. Therefore, since our BDM behaves as CDM 
away from the galactic center but the density has a core 
inner region, the proposed BDM profile is [i2[ 



Pbdm = pv,/[{r c/ts + r/rs){l + rjrsf 



(1) 



which contains three parameters: a typical scale length 
(ts) and density (po) of the halo, and a core radius ijc)- 
The BDM profile coincides with pnfw a-t large radius 
but has a core inner region, when the halo energy density 
Pbdm reaches the value pc = E"^ at r ~ Tc with 
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TABLE I: BDM parameters for different galaxies and mass 
model with rs^r^ in kpc, po in MQ/kpc^ and E^ in eV. 
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Pc = PBDM{r = Tc) ~ pors/2rc. 



(2) 



TABLE II: NFW parameters for different galaxies and 
mass models with same units as in given in Tab|l] 



The relevant parameters to determine the inner core 
structure of galaxies are the value of Ec = pj and r^- 
The parameters po and r^ depend on the initial condi- 
tions and formation of each galaxy. 

Analysis of LSB GalzLxies. Wc limit our sample to 
(early type and dwarf) THINGS galaxies with smooth, 
symmetric and cxtendcd-to-large-radii rotation curves. 
The set consists of seventeen low luminous LSB galax- 
ies. These observations represent the best available data 
to study the DM mass distribution and it has been used 
in works concerning the core versus cusp discrepancy con- 
troversy p. For technical details and systematic effects 
see [ij] . Our mass models include the three main compo- 
nents of a spiral galaxy: thin gaseous disk, Vgas , a thick 
stellar disk, T4, and a DM halo, Vhaio- We have not con- 
sidered a bulge because it is estimated to contain a small 
fraction of the total luminosity in our galaxy sample and 
therefore it has a small or null effect in the analysis. The 
gravitational potential of the galaxy is the sum of each 
mass component, thus the observed rotational velocity 
is Vobs'^ = V^^i^ + Vg^^, + Ti,V^. Avaflablc photometry, 
extracted from the SINGS images at the 3.6pm, shows 
that the stars in our sample of galaxies are distributed 
in a thin disk, with exponential central surface density 
profile J:{R) = Y.Qe'^/^'^. When Rd is the scale length of 
the disk and Sg is the central surface density with units 
[M0pc~^]. The rotation velocity of an exponential disk 
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FIG. 1: We show Tc vs Ec for Group A (circles ) and inner 
Group B (diamonds) galaxies, giving < Ec >~ 0.1 eV and core 
radius < Vc >~ 260 pc. Yellow (light) and blue (dark) represent 
(T and 2(T c.l. 



is given by the well know Freeman formula. The contri- 
bution of the gaseous disk is directly derived from the 
HI surface density distribution. The kinematics of stars 
brings a challenging problem in the analysis, mainly due 
to the uncertainty in the mass-to-light ratio (T^,), that is 
an additional constant, free parameter in the mass model. 
The stellar contribution to rotational curves is in many 
cases dominant close to the galactic center and this in- 
terferes with the determination of the core parameters. 
Some considerations has been made in order to reduce 
this uncertainty in the parameters [13| , but still the stel- 
lar contribution is not well known and depends on extinc- 
tion, star formation history, initial mass function (IMF), 
among other issues. We present a disk-halo decompo- 
sition using different assumptions for the stellar T^,: i) 
DM alone (Min.Disk); ii) DM and gas (Min.Disk-hgas); 
iii) Kroupa, considered as the minimal limit for the stel- 
lar disk; iv) diet-Salpeter, a stellar population synthesis 
model that yields a maximum stellar disk; and v) Free 
T^,, treating the mass-to-light ratio as an extra free pa- 
rameter in the model. 

We have analyzed the rotation curves of seventeen galax- 
ies using the above-mentioned five disk models and four 
different DM profiles (BDM, NFW, Burkert, ISO). Since 
we have a large number of data, and because we want to 
emphasize the main results, here we present the general 
conclusions and we show results for the Min.Disk scenario 
for BDM and NFW, only. However, a full comparison 
and a detailed analysis between all the different profiles 
and mass models can be found in [ij]. We also perform 
1 and 2 a likelihood contour plots for the BDM parame- 
ters, Tc and po for the different galaxies and mass models 
p^ . We have grouped the galaxies into three blocks ac- 
cording to the ratio rc/rg of the central fitted value for 
the min.disk scenario. The first group (G.A.) has a fitted 
Tc/rs < 1, the second group (G.B.) has rc/r's ~ 1, and fi- 
nally, the third group (G.C.) has rc/vg < 10^^. We shaU 
explain the physical interpretation of the results for each 
group in the following paragraphs. Since the energy Ec 
and Tc can take values (0, oo) the correct distribution for 
the sample that we considered is the log-normal distribu- 
tion, which is used to compute the statistics for the differ- 
ent BDM parameters. We obtain seven galaxies (DDO 
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TABLE III: BDM parameters and the slope a obtained 
from the fittings of the inner galactic region with the pi„ and 
pa profiles. The maximum distance is -Rmax(fcpc). 



154 (Al), NGC 2841 (A2), NGC 3031 (A3), NGC 3621 
(A4), NGC 4736 (A5), NGC 6946 (A6), and NGC 7793 
(A7)) (we call Group A) with fitted values rs > Vc ^ 
that are shown in Table ID For these galaxies we obtain 
average values for: the core Vc — 40 pc, Ts ^ 5 kpc, and 
the energy Ec ~ 0.11 eV for the min.disk analysis, that 
arc typical values for a galaxy, except for DDO 154 which 
has a larger core, Vc = 1.35 kpc. We notice that for all 
these galaxies the x^ in BDM is smaller than NFW's but 
the extra parameter (vc) in BDM makes the reduced Xred 
equivalent for both profiles. When comparing the BDM 
profile with the other two cored profiles, Burkert and ISO 
(Ref. [ij]) we conclude that BDM is still better. In par- 
ticular, Burkert and ISO profiles have difficulties when 
fitting a couple galaxies (NGC4736, NGC 3621) having 
xfed > 5 while BDM and NFW each have xled < 2. 
In the second case, Group B, (IC 2574 (Bl), NGC 2366 
(B2), NGC 2903 (B3), NGC 2976 (B4), NGC 3198 (B5), 
NGC 3521 (B6), and NGC 925 (B7)), we get a fitted value 
of Tc ~ r-g and BDM profile fits better than NFW. In fact, 
the NFW profile in some cases over-predicts the velocity 
in the inner parts of the galaxy or does fit with unphys- 
ical values for r, > 0(10^) or po < C'(lO3)M0/kpc3, see 
Table H] and ini The BDM profile gives r^ ~ r^ for all stel- 
lar mass models, having an average Tc value of '^6 and ~5 
kpc for Min.Disk and Min.Disk+Gas, respectively. Since 
Tc ~ Ts, implying that our BDM profile is of the type 
p (X {vc + r)~^ and hence having a strong core behavior. 
For this set of galaxies the Burkert and ISO profiles, as 
well as the BDM, fit much better than NFW, indicating 
that a core is needed. The BDM and Burkert profiles are 
slightly better than ISO when we consider Min.Disk and 
Min.Disk-|-Gas analysis. Otherwise, for Kroupa and diet- 
Salpeter stellar mass models the BDM and ISO profiles 
have equivalent Xred ^^'^ ^^ slightly better than Burk- 
ert. To study in more detail the central region of the 
galaxies, we performed an analysis of the inner galactic 
region [lj| . We define the inner region as the set of data 
for each galaxy with a constant slope value for the ob- 
served velocity, and for radius smaller than the fitted r^ 
in Table HI and before the slope of v starts to decrease 
[l3t . We computed the core radius, the central energy 
density, and the inner slope of the profile by examining 
the central region data with two different approaches: i) 
The first one by taking the limit r -C r^ in the BDM 
profile which gives pi„ = pors(rc -f r)^^, ii) The second 



is with the ansatz profile pa = Po'i^~°', where a is the 
slope a = — dlog[p]/log[r]. Both approaches are related 
by the equation < a ~ r/{rc + r) < 1, from which 
we obtained different values for a depending on the an- 
alyzed region, i.e. the slope in BDM takes the values 
a = (0, 1/2, 1) for r = {0,rc,rc < r < r,). The results 
arc in Table IIIII showing that in both approaches a core 
region is preferred and with a slope value in the interval 
< a < 0.52 for most of the galaxies and a maximum 
slope of a '^ 0.84, that implies that one is in the core 
dominant region, with r/r^ = «/(! — a), or at most in 
the transition region between Tc and Tg. We also obtain 
that in case i) the length of the core radius is reduced, 
the average Tc = 1.8 kpc with energy Ec ^ 0.06 eV for 
min.disk scenario, refer to Table lllll 
We only present here the conclusions of the last group C, 
composed of three galaxies (NGC2403, NGC50455 and 
NGC7331) [13 • These galaxies have fitted values such 
that Tc/rg < 10~^ for most stellar mass models. The 
reduced number of data close to galactic center prevents 
us from finding whether thesegalaxies actually have or 
not a core. We conclude in [13| that ISO and Burkert 
are poorly fitting profiles, meanwhile BDM has a central 
value for Tc ^' 0, so BDM reduces to NFW profile and 
both have equivalent fits. By applying inner analysis for 
these galaxies we obtain < Tc < 50 pc as a plausible 
interval that can fit the observations within the 2 a er- 
ror, although a more definitive conclusion would demand 
more data close to the center in these galaxies. 
We have also calculated fro the BDM profile the 1 and 
2 <T likelihood contour plots of Vc and po for the differ- 
ent galaxies and mass models, but due to lack of space 
here we refer to |13| . In general we found that when 
more mass components are included in the analysis the 
confidence level area becomes broader and in some cases 
the confidence levels increase up to an order of magni- 
tude for the Tc value [13 . Though Tc is different for each 
disk mass model there is an interval of values where Vc 
is consistent with all mass models within the 2 a errors. 
We also obtained that in the Free T^, scenario, the fitted 
value never gets a greater value than diet-Salpeter, which 
is consistent with [^. We show in Fig. [T]the 1 and 2 a 
contour plots of Ec vs r^ for the Min.Disk case and in Ta- 
bic IIVI shows the averages for the different galaxies and 
mass models and the 1 a standard deviation of Ec and 
Tc, for a detailed analysis see [l^- We found that when 
more mass components are taking into account, mainly 
when the contribution of the stellar disk has a dominant 
behavior close to the galactic center, similiar as in Ref. 
[18| , the evidence of the core fades away. We also found 
that even when re is different for each mass model, its 
value lays inside the confidence levels obtained from the 
other mass models for the same galaxy. 
Energy Phase Transition Ec- The values of the core 
radius re and the transition energy Ec = pj are shown 
in Table ID In Fig. [T] we plot Ec vs re obtained from the 
minimal disk analysis with its respectively contour errors. 
The circles and diamonds represent each galaxies from 



Group A and Group B, respectively. The energy of tran- 
sition between HDM and CDM takes place in Min.Disk 
mass model at energies and core radius up to 1 tr c.l. : 



Ec = 0.11 X io±0-46 eV, 



260 X 10=^'-^' pc. (3) 



The figures of Ec vs re for the other mass models 
(Min.Disk-|-Gas, Kroupa, and diet-Salpeter mass mod- 
els) are in [iJl. Notice that the dispersion on Ec is 
much smaller than that of re and this is consistent with 
our BDM model since Ec is a new fundamental scale 
for DM model. Even though we propose Ec as a new 
fundamental constant for DM it is important to stress 
that its value depends on the choice of BDM profile 
and on high resolution data for the rotation curves. We 
can also theoretically estimate the value of Ec and Ac 
which are related via pc = Ec = Ac n. We can then 
extract information from the underlying gauge group us- 
ing gauge group dynamics. Using Ac = A^ e~^^/^^' with 
Ac ~ £^c = 0.1 eV and gf = 1/2 the gauge coupling at 
unification or infiation scale A^ = 10^^ GeV. We obtain 
b = 8Tr/{gfLog[A,/Ac]) ci 1 and using that for a SUSY 
gauge group one has b = 3Nc — Nf we could have for 
example a gauge group SU{Ne — 3) with Nf = 8 funda- 
mental particles. The energy transition in Eq.([3|) is sim- 
ilar to the mean energy of a relativistic neutrino which 
is such that Y^m^ < 0.58 eV(95% CL) if one assumes 
3.04 neutrino species with degenerate mass eigenstates. 
Furthermore, an interesting connection could be further 
explored between the phase transition scale Ec of our 
BDM model, which also sets the mass of the BDM par- 
ticles, to the neutrino mass generation mechanism |19| . 
However, we would like to emphasize that our BDM are 
not neutrinos, since neutrinos are HDM, and a combina- 
tion of CDM plus neutrinos would have a cuspy NFW 
profile with CDM dominating in the inner region of the 
galaxies [3|. A main difference between neutrinos and 
BDM is that neutrinos were in thermal equilibrium at 
E ^ MeV with the standard model (SM) particles but 
BDM may not have been in thermal equilibrium with 
SM, as for example axion fields which are CDM but have 
a small mass, m < eV. 







BDM Statistics 








Energy Ec_ 


Core Tc 


Mass Models 


Ea„ 


Ea_ Be Sc+ 


O-Ee 


fc^ 


rc„ fc Ta^ 


(Trc 


Min.Disk 


0.14 


0.04 0.11 0.33 


0.46 


1.91 


0.01 0.26 5.41 


1.31 


M.Disk+gas 


0.09 


0.05 0.08 0.16 


0.27 


1.13 


0.12 0.63 3.26 


1.31 


Kroupa 


0.09 


0.02 0.06 0.21 


0.52 


6.92 


0.13 1.67 21.59 


1.11 


dict-Salpctcr 


0.03 


0.02 0.03 0.05 


0.24 


18.08 


3.22 12.27 46.68 


0.58 



TABLE IV: We show the statistics Ec(eV) and rc{kpc) for 
the BDM profile for the different mass models. In Columns (2) 
is the arithmetic mean energy, Ec„{eV), in columns (3-5) we 
present the median energy Ec and Ec± = EcW assuming a 
log-normal distribution. In columns (6-9) we present the same 
statistics for rc{kpc). 

Conclusions. To summarize, we have tested differ- 
ent DM profiles with several mass models using LSB 
THINGS galaxies. The study favors BDM (core pro- 
files) over the cuspy NFW, in accordance to general re- 



suits found in other works 0, m, 1201 ■ However different 
systematics may play an important role in the obser- 
vations such as noncircular motions, resolution of data 
and other issues p, Q. We presented here the analy- 
sis with Min.Disk for BDM and NFW profiles, and fur- 
ther we have performed the study using the five differ- 
ent disk mass components for the BDM, NFW, Burkert, 
and Pseudo-Isothermal profiles in |13j . We performed an 
analysis of the inner regions of Group B and C galaxies 
showing the ability to extract information from the rota- 
tion curves and the results are consistent with a core pro- 
file and we compute the core radius, central density and 
inner slope corresponding to the core region. In the cases 
where rc/rg < 10^^ the BDM is indistinguishable from 
NFW profile and much better than Burkert or ISO pro- 
files. The likelihood contour plots reveal how the central 
contribution of the stellar disk fades away the evidence 
of the core making Tc — >■ 0. The observational resolu- 
tion of the THINGS sample is of high quality, but one 
still needs data closer to the center of some galaxies, on 
scales smaller than 200 pc, in order to discern between 
cored or cuspy profiles. This is because stars pose a very 
challenging problem when testing the core-cusp problem, 
largely due to the uncertainty of the mass-to-light ratio 
and their dominant behavior close to the center of the 
galaxy. We computed the 1 and 2 a confidence levels for 
the BDM parameters Ec and Vc for the different galaxies 
and mass models. The contours plots show consistency in 
the fitted values of Ec and Vc for the different mass models 
for each galaxy. We found that the energy of transition 



between HDM and CDM takes places in Min.Disk mass 
models with an average energy Ec = Cll^Qgy eV and 
with a core radius re ^ 260 x lO*^'^^ pc, in accordance 
with Rcf. [21[. Confidence level contours show that the 
core radius depends on the number of disk mass compo- 
nents taking into account in the analysis. We found that 
Tc is highly constrained if the stellar disk has a dominant 
behavior close to the center of the galaxy. However, we 
notice that even though re is different for each disk mass 
model there is an interval of values where re is consistent 
for all mass models within the 2 a errors. The dispersion 
on Ec is much smaller than that of re and this is consis- 
tent with our BDM model since Ec is a new fundamental 
scale for DM model while re depends on r^ , po which de- 
pends on the formation and initial conditions for each 
galaxy. Given the richer structure that the BDM profile 
has, through its transition from CDM to HDM (when 
p ~ Pc), it allows for a physical explanation of the rota- 
tion curves of the different galaxies. 
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